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Abstract 


A search for the flavour-changing neutral-current decay t ^ qZ is presented. Data collected 
by the ATLAS detector during 2012 from proton-proton collisions at the Large Hadron Col¬ 
lider at a centre-of-mass energy of = 8 TeV, corresponding to an integrated luminosity 
of 20.3 fb”L are analysed. Top-quark pair-production events with one top quark decaying 
through the t ^ qZ {q = u,c) channel and the other through the dominant Standard Model 
mode t bW are considered as signal. Only the decays of the Z boson to charged leptons 
and leptonic W boson decays are used. No evidence for a signal is found and an observed 
(expected) upper limit on the t qZ branching ratio of 7 x 10”"^ (8 x 10”^) is set at the 95% 
confidence level. 
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1 Introduction 


The top quark is the heaviest elementary particle known, with a mass nit = 173.21 ± 0.51(stat.) + 
0.71(syst.) GeV [1]. Its lifetime is so short that, within the Standard Model (SM) of particle physics, 
it decays (almost exclusively to bW) before hadronisation occurs. These properties make it a particle well 
suited to test the predictions of the SM. In the SM, due to the GIM mechanism [2], flavour-changing neut¬ 
ral current (FCNC) decays such a&t^qZ are forbidden at tree level. They are allowed at one-loop level, 
but with a suppression factor of several orders of magnitude with respect to the dominant decay mode [3]. 
However, several SM extensions predict higher branching ratios (BRs) for the top-quark FCNC decays. 
Examples of such extensions are the quark-singlet model (QS) [4], the two-Higgs-doublet model with (FC 
2HDM) or without (2HDM) flavour conservation [5], the minimal supersymmetric model (MSSM) [6], 
supersymmetry with R-parity violation SUSY) [7] or models with warped extra dimensions (RS) [8]. 
For a review see Ref. [9]. The maximum values for the t ^ qZ BRs predicted by these models and by 
the SM are summarised in Table 1. Experimental limits on the ECNC t ^ qZ BR were established by 
experiments at the Earge Electron Positron Collider (EEP) [10-14], HERA [15], Tevatron [16, 17] and 
Earge Hadron Collider (EHC) [18, 19]. The most stringent limit, BR(t ^ qZ) < 5 x 10~^, is the one from 
the CMS Collaboration [19] using 25 fb“' of data collected at = 7 TeV and = S TeV. Previous 
ATEAS results obtained at = 7 TeV are also reported [18]. Eimits on other ECNC top-quark decay 
BRs (t —> qX, X = y,g, H) are reported in Refs. [10-14, 20-28]. 

This paper presents the ATEAS results from the search for the ECNC decay t ^ qZ in tt events produced 
at = 8 TeV, with one top quark decaying through the ECNC mode and the other through the SM 
dominant mode (f —> bW). Only the decays of the Z boson to charged leptons and leptonic W boson 
decays are considered. The final-state topology is thus characterised by the presence of three isolated 
charged leptons, at least two jets, and missing transverse momentum from the undetected neutrino. The 
paper is organised as follows. A brief description of the ATEAS detector is given in Section 2. The 
collected data samples and the simulations of signal and SM background processes are described in 
Section 3. Section 4 presents the object definitions, while the event analysis and kinematic reconstruction 
are explained in Section 5. Background evaluation and sources of systematic uncertainty are described in 
Sections 6 and 7. Results are presented in Section 8 and conclusions are drawn in Section 9. 


2 Detector and data samples 

The ATEAS experiment is a multi-purpose particle physics detector consisting of several sub-detector 
systems, which cover almost fully the solid angle' around the interaction point. It is composed of an in¬ 
ner tracking system close to the interaction point and immersed in a 2 T axial magnetic field produced by a 
thin superconducting solenoid, a lead/liquid-argon (EAr) electromagnetic calorimeter, an iron/scintillator- 
tile hadronic calorimeter, copper/EAr hadronic endcap calorimeter and a muon spectrometer with three 
superconducting magnets, each one with eight toroid coils. The forward region is covered by additional 
EAr calorimeters with copper and tungsten absorbers. The combination of all these systems provides 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the interaction point to the centre of the LHC ring, and the 
^-axis points upward. Cylindrical coordinates (r, (p) are used in the transverse plane, (p being the azimuthal angle around the 
beam pipe. The pseudorapidity is defined in terms of the polar angle d as rj = - lntan(6/2). The AR distance is defined as 
AR= + (A0F. 
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Table 1: Maximum allowed FCNC t ^ qZ BRs as predicted by several models [3-9] 


Model: 

SM 

QS 

2HDM 

LG 2HDM 

MSSM 

Ik SUSY 

RS 

BR(t ^ qZy. 

10-14 

10-^ 

10“® 

10-10 

10”^ 

10-® 

10-5 


charged-particle momentum measurements, together with efficient and precise lepton and photon identi¬ 
fication in the pseudorapidity range \t]\ < 2.5. Energy deposits over the full coverage of the calorimeters, 
It/I < 4.9 are used to reconstruct jets and missing transverse momentum (with magnitude A three- 

level trigger system is used to select interesting events. The Level-1 trigger is implemented in hardware 
and uses a subset of detector information to reduce the event rate to a design value of at most 75 kHz. 
This is followed by two software-based trigger levels which together reduce the event rate to less than 
1 kHz. A detailed description of the ATLAS detector is provided in Ref. [29] . 

In this paper the full 2012 dataset from proton-proton {pp) collisions ai ^[s - % TeV is used. The 
analysed events were recorded by single-electron or single-muon triggers and fulfil sfandard dafa-qualify 
requiremenfs. Triggers wifh differenl fransverse momenfum fhresholds are used fo increase fhe overall 
efficiency. The friggers using a low fransverse momenfum (px) threshold > 24 GeV) also have an 
isolation requirement. Efficiency losses at higher px values are recovered by higher threshold triggers 
(p^ > 60 GeV or > 36 GeV) without any isolation requirement. The integrated luminosity of the 
analysed data sample is 20.3 fb“^ 


3 Simulated samples 


In the SM, top quarks are produced at the LHC mainly in pairs, with a predicted tt cross section in pp col¬ 
lisions at a centre-of-mass energy of = 8 TeVof cr,f - 253]|]j5 pb for a top-quark mass of 172.5 GeV. 
The cross section has been calculated at next-to-next-to leading-order (NNLO) in QCD including resum¬ 
mation of next-to-next-to leading logarithmic (NNLL) soft gluon terms with Top-i-i- 2.0 [30-35]. The par- 
ton distribution function (PDL) and as uncertainties are calculated using the PDL4LHC prescription [36] 
with the MSTW 2008 68% CL NNLO [37, 38], CTIO NNLO [39, 40] and NNPDL 2.3 5f LLN [41] 
PDL sets and are added in quadrature to the renormalisation and factorisation scale uncertainties. The 
cross-section value for the NNLO-i-NNLL calculation is about 3% larger than the exact NNLO prediction 
implemented in HATHOR 1.5 [42]. 

The simulation of signal events is performed with PROTOS 2.2 [43, 44], which includes the effects of new 
physics at an energy scale A by adding dimension-six effective terms to the SM Lagrangian. The most 
general Ztu vertex that arises from the dimension-six operators can be parameterised including only 
and cr^’'^v terms [45] as: 


Lztu = 

where g is the electroweak coupling, cw is the cosine of the weak mixing angle, u and t are the quark 
spinors, is the Z boson field, Pl (Pr) is fhe lefl-handed (righf-handed) projection operafor, mz is fhe Z 
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boson mass and q'’ = p\ - is the outgoing Z boson momentum. The Ztc vertex can be parameterised 
in a similar fashion. This vertex involves a minimum of four anomalous couplings /<^p which 

are set to 0.01 each. It was checked that the coupling choice does not affect the kinematics of the event. 
No impact in the kinematics is seen by comparing the bWuZ and bWcZ processes and the latter is used 
as reference. Only decays of the W and Z bosons involving charged leptons are generated at the matrix- 
element level by PROTOS (Z ^ e^e~,p^p~, and W ev,pv, tv). The CTEQ6L1 [46] leading-order 
PDF is used. To account for higher-order contributions in the signal production, the events are reweighted 
to the measured tt differential cross section as a function of the transverse momentum of the tt system 
(l/cr)(dcr/dp") [47]. Hadronisation is handled by PYTHIA 6.426 [48] with the Perugia2011C [49] set 
of tuned parameters and r decays are processed with TAUOLA [50]. The top-quark mass is set to m, = 
172.5 GeV. Additional simulations with different parton shower parameterisations are used to estimate 
the systematic uncertainties on the amount of initial- and final-state radiation (ISR/FSR). 

Several SM processes have final-state topologies similar to the signal, with at least three prompt charged 
leptons, especially WZ, ZZ, ttV,^ ttH, ggH, VH, tZ and triboson {WWW, ZWW and ZZZ) production. 
Events with non-prompt leptons or in which at least one jet is misidentified as an isolated charged lepton 
(labelled as “fake leptons” throughout this paper) can also fulfil fhe even! selection requiremenfs. These 
evenfs, fypically Z-i-jels, Z + y, tt and single-fop, are esfimafed from a dafa-driven mefhod using a para- 
meferisafion of fhe frue- and fake-lepfon efficiencies. Samples of simulated evenfs of fhese backgrounds 
wifh fake lepfons are used fo cross-check fhe dafa-driven esfimafion. The Z-i-jefs simulafions include Z 
producfion in associafion wifh heavy-havour quarks. 

Table 2 summarises fhe informafion abouf fhe generafors, parfon shower and parfon disfribufion funclions 
used fo simulafe fhe differenl even! samples considered in fhe analysis. Diboson evenfs (VTZ and ZZ, 
where Z means Zjy*) produced using SHERPA confain up fo fhree addifional parfons and are selecfed fo 
have lepfons wifh pt > 5 GeV and mtt > 0.1 GeV for fhe Zjy*. The addifional WZ and ZZ samples are 
used for comparison. The WZ ALPGEN samples are simulafed wifh up fo five addifional parfons from fhe 
mafrix elemenf. The ZZ HERWIG [51] samples are selecfed fo have one lepfon wifh pj > 10 GeV and 
\q\ < 2.8. The simulations of ttZ, ttW{W), tZ and tribosons include events with up to two extra partons in 
the final state. The simulated samples used to cross-check the data-driven estimation of background with 
fake leptons are also listed in Table 2. 

Detailed and fast simulations of the detector and trigger are performed with standard ATEAS software 
using GEANT4 [52, 53] and ATLFASTII [53], respectively. The same offline reconstruction methods used 
on data are applied to the simulated samples. 


4 Object reconstruction 


The primary physics objects considered in this search are electrons, muons, jets, and fj-tagged jets. 
Tau leptons are not explicitly reconstructed, although the t decay products are reconstructed as electrons, 
muons or jets and as an additional contribution to the missing transverse momentum. 

Electron candidates are reconstructed [63] from energy deposits (clusters) in the electromagnetic calor¬ 
imeter, which are then matched to reconstructed charged-particle tracks in the inner detector. The can¬ 
didates are required to have a transverse energy Et greater than 15 GeV and a pseudorapidity of the 
calorimeter cluster associated with the electron candidate |?7ciusterl < 2.47. Candidates in the transition 

^ V stands for W and Z bosons. The simulation of tiW events also includes a possible second W boson. 
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Table 2; Generators, parton shower, parton distribution functions and parameter tune for hadronisation used to 
produce simulated samples used in this analysis. 


Sample 

Generator 

Parton shower 

PDF 

Tune 

Samples with at least three prompt leptons: 




bWqZ 

PROTOS 2.2 [43] 

PYTHIA 6.426 [48] 

CTEQ6L1 [46] 

Perugia2011C [49] 

WZ 

SHERPA 1.4.3 [54] 

SHERPA 1.4.3 

CTIO [39] 

— 

WZ 

ALPGEN 2.14 [55] 

HERWIG 6.520.2 [51] 

CTEQ6L1 

AUET2 [56] 

zz 

SHERPA 1.4.3 

SHERPA 1.4.3 

CTIO 

— 

zz 

HERWIG 6.5 

HERWIG 6.5 

CTEQ6L1 

AUET2 

ttV, tZ, tribosons 

MADGRAPH5 1.3.33 [57] 

PYTHIA 6.426 

CTEQ6L1 

AUET2B 

tlH, WH, ZH 

PYTHIA 8.163 [58] 

PYTHIA 8.163 

CTEQ6L1 

AU2 [59] 

ggfi 

POWHEG 2 [60] 

PYTHIA 8.163 

CTIO 

AU2 

Other samples: 






SHERPA 1.4.3 

SHERPA 1.4.3 

CTIO 

— 

Z+jets (30 GeV < < 1 TeV) 

ALPGEN 2.14 

PYTHIA 6.426 

CTEQ6L1 

Perugia2011C 

Z-f-jets (10 GeV < mce < 60 GeV) 

ALPGEN 2.14 

HERWIG 6.520.2 

CTEQ6L1 

AUET2 

Zy 

SHERPA 1.4.1 

SHERPA 1.4.1 

CTIO 

— 

tt^ bWbW 

POWHEG 2 

PYTHIA 6.426 

CTEQ6L1 

Perugia2011C 

Single top (s, Wt channel) 

MC@NL0 4.03 [61] 

HERWIG 6.520.2 

CTIO 

AUET2 

Single top (t channel) 

AcerMC 3.8 [62] 

PYTHIA 6.426 

CTEQ6L1 

AUET2B 


region between the barrel and endcap calorimeters with 1.37 < Itydusterl < 1-52 are excluded. Electron 
candidates in this analysis must satisfy tight quality requirements on the electromagnetic cluster and as¬ 
sociated track which provide discrimination between isolated electrons and jets. In order to suppress 
multi-jet backgrounds, it is also required that there is little activity in the space surrounding the electron. 
Two isolation variables are employed: the energy deposited around the electron in the calorimeter in a 
cone of size AT? = 0.2 and the scalar sum of the pj of the tracks within cone of size AT? = 0.3 around the 
electron. Cuts on these two quantities are used to select isolated electrons; the adopted cuts yield a 90% 
identification efficiency in Z boson decays to events from the full 2012 dataset. Additionally, the 
longitudinal impact parameter |zol of the electron track with respect to the selected primary vertex of the 
event is required to be less than 2 mm. The closest jet if separated by AT? < 0.2 from the selected elec¬ 
tron is removed from the event. The electron candidate is discarded if an additional selected jet is found 
with AT? < 0.4. A looser electron selection, used for the estimation of backgrounds with fake leptons, is 
defined by removing fhe isolation requiremenfs. 

The muon candidafe reconstrucfion [64] is performed by finding, combining and tiffing frack segmenfs 
in fhe layers of fhe muon chambers, sfarfing from fhe oufermosf layer. The idenfitied muons are fhen 
mafched wifh fracks reconsfrucfed in fhe inner defector. The candidafes are retiffed using fhe complefe 
frack informafion from bofh defecfor sysfems, and are required fo safisfy pj > 15 GeV, \t]\ < 2.5 and 
fo be separafed by AT? > 0.4 from any selecfed jef. The hif paffem in fhe inner defecfor is required fo 
be consisfenf wifh a well-reconsfrucfed frack and fhe |zol of fhe muon track is required to be less than 
2 mm. Additionally, the sum of the momenta of tracks inside a cone around the muon candidate, with 
variable size such that it is smaller for higher muon p^ [65], must be less than 5% of the muon energy. For 
the estimation of backgrounds with fake leptons, a looser selection is applied by removing the isolation 
requirement. 

Jets are reconstructed [66] from topological clusters of neighbouring calorimeter cells with significant 


5 










energy deposits using the anti-^j algorithm [67] with a radius parameter R = 0.4. Prior to jet finding, a 
local calibration scheme is applied to correct the topological cluster energies for the non-compensating 
response of the calorimeter, dead material and energy leakage. The corrections are obtained from simu¬ 
lations of charged and neutral particles. These jets are then calibrated to the hadronic energy scale using 
Pt- and ? 7 -dependent correction factors. Dedicated requirements are applied to remove the negligible 
fraction of events (less than 0.01%) where a jet is incorrectly reconstructed from a few noisy calorimeter 
cells [68]. The jets used in the analysis are required to have px > 25 GeV and |p| < 2.5. To reduce the 
number of selected jets that originate from secondary pp interactions, for jets with px < 50 GeV and 
\t]\ < 2.4, the scalar sum of the px of tracks matched to a jet and originating from the primary vertex must 
be at least 50% of the scalar sum of the px of all tracks matched to the jet. 

Jets containing f^-hadrons are identified (‘fj-fagged’) [69] using an algorifhm based on mulfivariafe fech- 
niques. If combines informafion from fhe impacf paramefers of displaced fracks and from topological 
properties of secondary and tertiary decay vertices reconstructed within the jet. It is determined with 
simulated tt events that, for the chosen working point, the tagging efficiency for f^-jets with px > 20 GeV 
is 70%, while the rejection factors for light-quark or gluon jets (light jets), charm jets and t leptons are 
137, 5 and 13, respectively. 

The measurement of is based [70] on the energy deposits in the calorimeter with \r]\ < 4.9. The en¬ 
ergy deposits associated with reconstructed jets and electrons are calibrated accordingly. Energy deposits 
not associated with a reconstructed object are calibrated according to their energy sharing between the 
electromagnetic and hadronic calorimeters. The momentum associated with each reconstructed muon, 
estimated using the momentum measurement of its reconstructed track, is taken into account in the cal¬ 
culation ofE™"^ 

5 Event selection and kinematics 

At least one of the selected leptons must be matched, with AR < 0.15 to the appropriate trigger object 
and have px > 25 GeV. The trigger efficiencies for the leptons are approximately 93% for electrons, 70% 
for muons with Ipl < 1.05 and 86% for muons with 1.05 < Ipl < 2.4 [71, 72]. The events are required to 
have at least one primary vertex with more than four associated tracks, each with px > 400 MeV. The 
primary vertex is chosen as the one with the highest X Pj over all associated tracks. Leptons from cosmic 
rays are rejected by removing muon pairs with large, oppositely signed transverse impact parameters 
(|(fol > 0.5 mm) and consistent with being back-to-back in the r - <p plane. Events with noise bursts and 
readout errors in the LAr calorimeter are also rejected. Exactly three isolated leptons associated with the 
same vertex are required. The three leptons must have Ipl <2.5 and px > 15 GeV. Two of the leptons 
are required to have the same flavour, opposife charge and a reconsfrucfed mass wifhin 15 GeV of fhe Z 
boson mass (niz) [1]. If more fhan one compatible leplon-pair is found, fhe one wifh fhe reconsfrucfed 
mass closesf to mz is chosen as fhe Z boson candidafe. According fo fhe signal fopology, fhe evenfs 
are fhen required fo have > 20 GeV and fwo jefs, allhough an addilional Ihird jel from inilial- or 
final-slale radiafion is allowed. All jefs are required fo have px > 35 GeV and \r]\ < 2.5. One or fwo of 
fhe jefs musf be b-tagged. Only one b-tagged jef is expecfed in fhe signal evenfs, neverfheless a second 
one can arise from a misidenfified c-jef associaled wifh fhe ECNC decay of fhe top quark. Allowing for 
fhe additional (^-tagged jel increases fhe signal efficiency wilhout compromising fhe signal-fo-background 
rafio. 
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Applying energy-momentum conservation, the kinematics of the top quarks can be reconstructed from 
the corresponding decay particles. Since the neutrino from the semileptonic decay of the top quark 
(t ^ bW —> bbv) is undetected, its four-momentum must be estimated. This can be done by assuming that 
the lepton not previously assigned to the Z boson and the Z^-tagged jet (labelled b-]&t) originate from the 
W boson and SM top-quark decays, respectively, and that is the neutrino’s transverse momentum. 
The longitudinal component of the neutrino’s momentum (pp is then determined by minimising, without 
constraints, the following expression: 


A 


2 




( 2 ) 


where and are the reconstructed masses of the qZ, bW and bv systems, respectively. 

The central value for the masses and the widths of the top quarks and W boson are taken from recon¬ 
structed simulated signal events. This is done by matching the true particles in the simulated events to 
the reconstructed ones, setting the longitudinal momentum of the neutrino to the of the true simulated 
neutrino and then performing Bukin fits^ [73] to the masses of the matched reconstructed top quarks 
and W boson. The values are = 173 GeV, = 10 GeV, - 168 GeV, - 23 GeV, 

mvr - 82 GeV and crw - 15 GeV. 

For each jet combination, where any jet can be assigned to ja, while ji, must correspond to a (j-tagged 
jet, the minimisation gives the most probable value for pp From all combinations, the one with the 
minimum is chosen, along with the corresponding p\ value. The jet from the top-quark FCNC decay 
is referred to as the light-quark {q) jet. The fractions of correct assignments between the reconstructed 
top quarks and the true simulated particles (evaluated as a match within a cone of size A/? = 0.4) are 
ffPCNc - 79.9% and = 56.3%. Figure 1 shows the px of the third lepton as well as the and;^^ 
distributions at this level of the analysis. 

The selection of the signal region is concluded with the requirement of x^ < 6, which optimises the 
sensitivity discussed in Section 8. 


6 Background estimates 


Three control regions are defined to check the agreement between data and simulated samples of the ZZ, 
VTZ and ttZ backgrounds. No scaling factors are derived from these control regions, however they are 
used to estimate the background modelling uncertainties described in Section 7. The tZ contribution to 
the total background is expected to be smaller than the one from ttZ events [74]. Due to the similarity 
between the final sfafes of tZ and signal evenfs, fhere are large signal confribufions fo possible tZ confrol 
regions. For fhese reasons no confrol region is defined for fhe tZ background. 

The ZZ confrol region is defined by requiring fwo pairs of leptons wifh fhe same flavour, opposife charge 
and a reconsfrucfed mass wifhin 15 GeV of fhe Z boson mass. The expecfed and observed yields are 
shown in Table 3 and fhe SHERPA sample is chosen as reference. 

To sfudy fhe fkZ background fhe following confrol region is defined. Evenfs are required fo have fhree 
leptons, fwo of fhem wifh fhe same flavour, opposite charge and a reconsfrucfed mass wifhin 15 GeV 
of fhe Z boson mass. Additional requiremenfs include fhe presence of af leasf one jef in fhe even! wifh 

^ This is a piecewise function with a Gaussian function in the centre and two asymmetric tails. 
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Figure 1: Expected (filled histogram) and observed (points with error bars) distributions in the signal region before 
the cut is applied for (a) px of the third lepton, (b) E™** and (c) x'^ of the kinematics reconstruction. For 
comparison, distributions for the FCNC tt —> bWqZ signal (dashed line), normalised to the observed 95% CL limit 
reported in this paper, are also shown. Background statistical uncertainties associated with the number of events in 
the samples are represented by the hatched areas. 


Table 3: Event yields in the ZZ control region for all significant sources of background. The ZZ SHERPA sample is 
taken as reference for the total background estimation. The first uncertainty is the statistical one associated with the 
number of events in the simulated samples, the second uncertainty is systematic and is described in Section 7. The 
entry labelled “other backgrounds” includes all the remaining backgrounds described in Section 3 and in Table 2. 
The signal efficiency is also shown. 


Sample 

Yields 

ZZ (SHERPA) 

87 ± 4 ± 5 

ZZ (HERWIG) 

85±4±5 

Other backgrounds 

0.48 ± 0.05 ± 0.08 

Total background 

88±4±5 

Data 

95 

Signal efficiency [xlO“**] 

5.6 ±4.3 ±0.1 










































Table 4: Event yields in the WZ control region for all significant sources of background. The VTZ SHERPA sample 
is taken as reference for the total background estimation. The hrst uncertainty is the statistical one associated with 
the number of events in the samples, the second uncertainty is systematic and is described in Section 7. The entry 
labelled “other backgrounds” includes all the remaining backgrounds described in Section 3 and in Table 2. The 
signal efficiency is also shown. 


Sample 

Yields 

IFZ (SHERPA) 

333 + 5 ± 17 

IVZ (ALPGEN) 

393 ± 6 ± 19 

ZZ 

35 ± 3 ± 6 

Fake leptons 

15 ±3 ±5 

Other backgrounds 

9.5 ± 0.3 ± 2.4 

Total background 

392 + 7 ± 19 

Data 

405 

Signal efficiency [xl0“‘*] 

9.8 ±0.1 ± 1.0 


Table 5: Event yields in the tiZ control region for all significant sources of background. The first uncertainty is 
the statistical one associated with the number of events in the samples, the second uncertainty is systematic and is 
described in Section 7. The entry labelled “other backgrounds” includes all the remaining backgrounds described 
in Section 3 and in Table 2. The signal efficiency is also shown. 


Sample 

Yields 

ttV 

8.3 ± 0.2 ± 2.7 

tZ 

2.0 ±0.1 ± 1.0 

wz 

1.8 ±0.3 ±0.4 

Other backgrounds 

1.8 ±0.4 ±0.4 

Total background 

13.9 ±0.6 ±3.0 

Data 

12 

Signal efficiency [xl0~‘^] 

3.9 ±0.1 ±0.6 


P'Y > 35 GeV, no Z^-tagged jets with pj > 35 GeV and a W boson transverse mass, built with the 
residual lepton and greater than 50 GeV. Table 4 shows the expected and observed yields in this 

control region. The best estimation comes from the SHERPA prediction, which is chosen as the reference 
sample. 

The tiZ control region is defined by requiring at least three leptons, two of them with the same flavour, 
opposite charge and a reconstructed mass within 15 GeV of the Z boson mass. Furthermore the events are 
required to have at least two jets with pj > 25 GeV and at least two Z>-tagged jets if there are three leptons 
in the event, or at least one ^-tagged jet if there are four or more leptons in the event. Since the signal 
contribution for events with three leptons and two Z^-tagged jets is small, the overlap between signal and 
background regions is not removed, increasing the ttZ sensitivity in this control region. Table 5 shows 
the yields in this control region, and the background yields agree very well with the data within the given 
uncertainty. 

Backgrounds from events which contain at least one fake lepton are estimated from data using the matrix 
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Table 6: Event yields in the fake-lepton control region for all significant sources of background. The first uncertainty 
is the statistical one associated with the number of events in the samples, the second uncertainty is systematic and 
is described in Section 7. The entry labelled “other backgrounds” includes all the remaining backgrounds described 
in Section 3 and in Table 2. The signal efficiency is also shown. 


Sample 

Yields 

Fake lepfons 

7 ± 1 ±4 

WZ 

2.7 ± 0.4 ± 0.7 

ZZ 

1.7 ±0.6 ±0.8 

Ofher backgrounds 

1.7 ±0.1 ±0.6 

Tofal background 

13 ± 1 ±4 

Dafa 

17 

Signal efficiency [xl0“^] 

1.77 ±0.06 ±0.20 


method [75]. This is based on the measurement of the efficiencies of real and fake loose leptons to pass the 
nominal selection, and eF^ and on the selection of two orthogonal sets of events in the signal region. For 
the first of these sets, the nominal requirements are used for the leptonic selection, while for the second 
one, only the leptons which satisfy the looser selection (as described in Section 4) but without meeting 
the nominal requirements are considered. For the single-lepton case, the number of events with one fake 
nominal lepton is = (ep/fR - ep) ((fR - 1 )A^t + where Nj (Ni) represents the number of 

selected events in the first (second) set defined above. The mefhod is exfrapolafed fo fhe fhree-lepfon 
fopology, wifh a 8 x 8 mafrix fhaf is inverted using a numerical mefhod fo obfain fhe number of evenfs 
wifh af leasf one fake lepfon. The efficiencies for real and fake lepfons are esfimafed as a funclion of 
fhe lepfon fransverse momenfum by a fif of fhe mafrix mefhod resulfs fo fwo dedicafed enriched samples 
of real and fake lepfons: a sample of Z —> ,{ = e,iu and a same-sign dilepfon sample (excluding 

same-flavour evenfs wifh a reconsfrucfed mass compafible wifh a Z boson). In bofh samples, in order fo 
improve fhe modelling of fake lepfons originafing from heavy-flavour decays, only evenfs wifh af leasf 
one additional (j-fagged jef are considered. The efficiency fR ranges from 0.74 fo 0.88 (0.80 fo 0.99) and 
ep from 0.010 fo 0.13 (0.035 fo 0.18) for elecfrons (muons). The relevanf uncerfainfies are calculated from 
fhe discrepancy befween predicted and observed number of evenfs in fhe confrol region defailed below. 

A confrol region fo fesf fhe performance of fhe fake-lepfon estimation mefhod and derive ifs uncerfainfy 
is defined. If requires fhree lepfons wifh pj < 50 GeV (fhe fhird one wifh pj < 30 GeV), fwo of fhem 
having fhe same flavour, opposife charge and a reconsfrucfed mass wifhin 15 GeV of fhe Z boson mass, 
af leasf one f7-fagged jef wifh pj > 35 GeV and < 40 GeV. As for fhe ttZ confrol region, fhere is 
a small overlap wifh fhe signal region, which is nof removed in order fo increase fhe sensifivify fo fhe 
fake-lepfon backgrounds. The yields are shown in Table 6 and agree well befween dafa and expecfafion. 
As a validafion of fhe mafrix mefhod, fhe background in which exacfly one of fhe lepfons is a fake lepfon 
is also evaluated using simulated samples. The resulfs for fhe signal region and differenf confrol regions 
are consisfenf befween fhe fwo mefhods wifhin fhe esfimafed uncerfainfies. 

Figure 2 shows fhe pj of fhe leading lepfon for fhe ZZ, IFZ and ttZ confrol regions, and fhe reconsfrucfed 
mass of fhe fwo lepfons wifh fhe same flavour and opposife charge for fhe fake-lepfon confrol region. 
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Figure 2: Expected (filled histogram) and observed (points with error bars) distributions for the pj of the leading 
lepton in the (a) ZZ, (b) VFZ and (c) ffZ control regions and (d) reconstructed mass of the two leptons with the 
same flavour and opposite charge in the fake-lepton control region. For comparison, distributions for the FCNC 
tl bWqZ signal (dashed line), scaled to 10"^ or 10 times the observed 95% CL limit, are also shown. Background 
statistical uncertainties associated with the number of events in the samples are represented by the hatched areas. 
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7 Systematic uncertainties 


The effect of each source of systematic uncertainty is studied hy independently varying the corresponding 
central value and propagating this through the full analysis chain. The relative impact of each type of 
systematic uncertainty on the total background and signal is summarised in Table 7. 

The main uncertainty on the backgrounds comes from their modelling, which has the following two 
contributions. The level of agreement with data of the reference samples is assessed from the combination 
of the Poisson uncertainty on the available amount of data with the statistical uncertainty of the expected 
background yields in the dedicated control regions described in Section 6. The uncertainties are estimated 
to be 6.3%, 12%, 42% and 62%, for the fTZ, ZZ, ttZ and fake-lepton backgrounds, respectively. The other 
contribution comes from the uncertainty on the theoretical prediction in the signal region and is estimated 
using the alternative VkZ and ZZ simulated samples. The corresponding uncertainties are 17% and 100%, 
respectively. Similarly, for ttZ, tZ and Higgs samples, conservative values of 30% [76, 77], 50% [74] and 
15% [78] respectively, are used, in order to account for the theoretical uncertainties. The combination of 
all these uncertainties gives a 17% effect on the total background estimation. 

The theoretical uncertainties of the signal modelling, as described in Section 3, namely production cross 
section and ISR/FSR modelling, are found to be 5.5%. 

For both the estimated signal and background event yields, experimental uncertainties resulting from 
detector effects are considered. The lepton reconstruction, identification and trigger efficiencies, as well 
as lepton momentum scales and resolutions [63, 79, 80] are considered. The overall effect on the total 
background yield and the signal efficiency is estimated to be 4.7% and 2.9% respectively. The effect of the 
jet energy scale and resolution [66, 81] uncertainties are evaluated as 7.7% and 4.9% for the background 
and signal, respectively. The l^-tagging performance component, which includes the uncertainty of the b-, 
C-, mistagged- and r-Jet scale factors (the r and charm uncertainties are highly correlated and evaluated 
as such) is evaluated by varying the rj-, pj- and flavour-dependent scale factors applied to each jet in the 
simulated samples. It is estimated to be 3.9% for the total background and 7.2% for the signal efficiency. 
The scale uncertainty [70] is found to vary the total background yield and the signal efficiency by 
3.2% and 1.5%, respectively. All these detector systematic uncertainties are treated as fully correlated 
between signal and background. 

The uncertainty related to the integrated luminosity for the dataset used in this analysis is 2.8%. It is 
derived following the methodology described in Ref. [82]. It only affects the estimations obtained from 
simulated samples, therefore its impact on the total background yield estimation is 2.4%. 


8 Results 

Table 8 shows the expected number of background events, number of selected data events and signal 
efficiency after the final event selection described in Section 5. Figure 3 shows the reconstructed masses 
of the top quarks and Z boson after the final selection. Good agreement between data and background 
yields is observed at all stages of the analysis. No evidence for the t ^ qZ decay is found and a 95% CL 
upper limit on the number of signal events is derived using the modified frequentist likelihood method [83, 
84]. 
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Table 7: Summary of the impact of each type of uncertainty on the total background and signal yields. The values 
are shown as the relative variations from the nominal values. The statistical uncertainty associated with the number 
of events in the simulated samples is also shown. 


Source 

Background [%] 

Signal [%] 

Background modelling 

17 

— 

Signal modelling 

— 

5.5 

Leptons 

A.l 

2.9 

Jets 

7.7 

4.9 

f7-tagging 

3.9 

7.2 

^miss 

3.2 

1.5 

Luminosity 

2.4 

2.8 

Statistical 

8.1 

1.5 


Table 8: Expected number of background events, number of selected data events and signal efficiency (normalised 
to all decays of the W and Z bosons), after the final selection. The first uncertainty is the statistical one associated 
with the number of events in the samples, the second uncertainty is systematic and is described in Section 7. The 
entry labelled “other backgrounds” includes all the remaining backgrounds described in Section 3 and in Table 2. 


Sample 

Yields 

WZ 

1.3 ±0.2 ±0.6 

tW 

1.5 ±0.1 ±0.5 

tz 

1.0 ±0.1 ±0.5 

Fake leptons 

0.7 ± 0.3 ± 0.4 

Other backgrounds 

0.2 ±0.1 ±0.1 

Total background 

4.7 ± 0.4 ± 1.0 

Data 

3 

Signal efficiency [xl0“"^] 

7.8 ±0.1 ±0.8 


The test-statistic X^i, which compares the number of observed data events with background and signal 
expectations, is defined as: 

Xrf = nln(l + 0 (3) 

where n, s and b are the numbers of data, expected background and signal events, respectively. The 
statistical test is then compared to 10^ pseudo-experiments for the hypotheses of signal plus background 
(Xs+b) and background-only (Xb), which are obtained by replacing n with the corresponding number 
of events produced by each pseudo-experiment. The statistical fluctuations of the pseudo-experiments 
are implemented assuming that the number of events follows a Poisson distribution. All statistical and 
systematic uncertainties on the expected backgrounds and signal efficiencies, as described in Section 7, 
are taken into account and implemented assuming Gaussian distributions. 
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Figure 3: Expected (filled histogram) and observed (points with error bars) distributions in the signal region after 
the final selection is applied for the reconstructed masses of the (a) top quark from the FCNC decay, (b) top quark 
from the SM decay and (c) Z boson. For comparison, distributions for the FCNC tl —> bWqZ signal (dashed line), 
normalised to the observed 95% CL limit, are also shown. Background statistical uncertainties are represented by 
the hatched areas. 


The CL for a given signal hypothesis 5 is defined as [83]: 


1 - CL = 


jf" Ps+biX)dX 


(4) 


where Ps+t and Pt are the probability density functions obtained from the pseudo-experiments for the 
Xs+b and Xb values, respectively, and are functions of s and b. The limit on the number of signal events is 
determined by finding the value of s corresponding to a CL of 95%. The expected limit is computed by 
replacing X^ with the median of the statistical test for the background hypothesis {Xb). 

The limits on the number of signal events are converted into upper limits on the t ^ qZ branching fraction 
using the NNLO-i-NNLL calculation, and uncertainty, for the tt cross section, and constraining BR(t ^ 
bW) = 1 - BR(t ^ qZ). Table 9 shows the observed limit on BR(t —> qZ) together with the expected 
limit and corresponding ±lcr bounds. These values are calculated using the reference tt —> bWcZ sample, 
since it gives a more conservative result than the tt bWuZ sample. The smaller ^-tagged jet multiplicity 
of the ti bWuZ signal sample leads to an improvement of 4% in the limit. 

Figure 4 compares the 95% CL observed limit found in this analysis with the results from other FCNC 
searches performed by the HI, ZEUS, LEP (combined results of the AEEPH, DEEPHI, E3 and OPAE 
collaborations), CDE, D0 and CMS collaborations. The results presented in this paper are consistent with 
the ones from the CMS Collaboration. 
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Table 9: Observed and expected 95% CL limits on the FCNC top-quark decay BRs. The expected central value is 
shown together with the +lcr bands, which includes the contribution from the statistical and systematic uncertain¬ 
ties. 


Observed 

7 X 10-^ 

(-Icr) 

6 X 10-4 

Expected 

8 X 10-"^ 

(-1-1 cr) 

12 X 10-4 



(a) (b) 


Figure 4; The current 95% CL observed limits on the (a) BR(f ^ qy) vs BR(f —> qZ) and (b) BR(f —> qH) vs 
BR(f —> qZ) planes are shown. The full lines represent the results from the ATLAS [25], CDF [16, 21] CMS [19, 
27], D0 [17], HI [20], LEP (combined results of the ALEPH, DELPHI, L3 and OPAL collaborations) [10-14] and 
ZEUS [15] collaborations. The ATLAS lines correspond to the limit on BR(f —> qZ) set in this paper. 
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9 Conclusions 


A search for the FCNC top-quark decay t ^ qZ in events with three leptons has been performed using 
LHC data collected by the ATLAS experiment at a centre-of-mass energy of = 8 TeV and correspond¬ 
ing to an integrated luminosity of 20.3 fb“^ recorded in 2012. No evidence for signal events is found and 
a 95% CL limit for the t ^ qZ branching fraction is established at BR(t —> qZ) < 7 x 10~^, in agreement 
with the expected limit of BR(t ^ qZ) < 8 x 10““^. 
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